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I. 

We consider the Weyl-Dirac theory within the framework of the weak field approximation 
and show that the resulting gravitational potential differs from that of Newtonian by a repulsive 
correction term increasing with distance. The scale of the correction term appears to be determined 
by the time variation rate of the gravitational coupling. It is shown that if the time variation rate 
of gravitational coupling is adopted from observational bounds, the theory can explain the rotation 
£■ — ' curves of typical spiral galaxies without resorting to dark matter. To check the consistency of our 

theoretical model with observation we use Likelihood analysis to find the best-fit values for the 
free parameters. The mean value for the most important free parameter, (3 x 10 14 (l/yr), using 

+ V.O Tnn-Hat onrl Honoaian nrinrQ are R 00 +2.44+6.18 j r + l .22+2.90 



the Top-Hat and Gaussian priors are 6.38_ 3 ; 4 gl 6 ; 71 and 5.72_ 1 ; 18 I 2 ; 6g , respectively. Although the 
interval for which (3 is defined is wide, our results show that the goodness of the fit is, by and 
large, not sensitive to this quantity. The intergalactic effects and gravitational lensing of clusters 



of galaxies are estimated and seem to be consistent with observational data. 

\ 



^jrj; 1 Introduction 

It has long been known that standard gravitational theories cannot correctly predict the dynamics 
of large astronomical systems. The rotation curves of spiral galaxies and gravitational lensing are 
among the well known examples [1, 2]. The rotation velocities usually tend to a constant or slightly 
rising values as a function of distance from the center of galaxy. This is in sharp contrast to the 
inverse square force law which implies a decline in velocity. The gravitational lensing of clusters of 
galaxies is the other observation that appears to be in conflict with the standard gravitational theories 
on extra galactic scales [2]. These observations are usually explained by postulating the existence of 
dark matter [3, 2]. Since there is no established observation of dark matter, the belief that standard 
gravitational theories should be modified on large length scales has been gaining momentum in recent 
years. Many attempts have emerged in this regard recently [4]-[14]. 

The purpose of the present paper is to undertake an investigation on the predictions of the 
dynamics of large astronomical systems in the Weyl-Dirac theory [15] which, in effect, is the Weyl 
theory [15] modified by Dirac later on. Dirac made use of the theory to provide a framework to 
explain his large number hypothesis. This modified theory can be considered as a scalar-tensor 
theory of gravity and leads to a time varying gravitational coupling G in accordance with the Dirac 
hypothesis. In comparison with Weyl theory, this modified theory is simpler and in agrement with the 
general theory of relativity. As we shall show below, its implications on the dynamics of astronomical 
objects is significant and therefore well worth studying. 
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The organization of this paper is as follows: in the next section we briefly discuss the ingredients 
of the Weyl geometry followed by the introduction of the Weyl-Dirac model in section 3. In section 
4 we present the spherically symmetric solutions of the model for a point-like mass source [16]. In 
section 5 we generalize the solutions to galaxies considered as spherically symmetric extended mass 
sources. The consistency of our model with observational rotation curves of typical spiral galaxies 
using Likelihood statistics will be checked and the best-fit values for free parameters of the model are 
determined in section 6. Sections 7 and 8 are devoted to the intergalactic effects and the gravitational 
lensing. Concluding remarks are presented in section 9. 

In what follows we shall use units in which h = c = G = 1. However, in order to interpret 
physical results we turn to conventional units from time to time. The signature of the metric is taken 
as (-,+,+,+) and the convention for the sign of the curvature tensor is that of Misner, Thorne and 
Wheeler [17]. For simplicity, the scalar field indices are taken as their covariant derivatives. 

2 Weyl geometry 

Weyl geometry is a natural generalization of the Riemannian geometry [15]. In Weyl geometry the 
length of a vector is assumed to change as well as its direction in a parallel displacement. This means 
that if a vector has length I at a point with coordinates then 

SI = (k^)l , (1) 

after a parallel displacement Sx^. The covariant vectors k^ are considered as the field quantities and 
may be called the Weyl meson fields. By equation (1), the total displacement around a small closed 
loop is 

SI = IF^SS^ , (2) 

where 

F[w = kv,fi ^1,1/ , (3) 

is called length curvature [18] and SS^ U is the element of the area enclosed by the loop. It follows 
from definition (1) that the comparison of elements of length at two different points which are not 
separated by an infinitesimal distance can be made only with respect to a path joining the points and 
obviously, different paths lead to different results. The comparison then is possible if one arbitrarily 
defines the standards of length at each space-time point. Alternatively, the change in length I can be 
parameterized by means of an arbitrary function 0,(x^) such that 

I' = n(x^)l , (4) 

and 

^ = ^ + (ln[^)])^. (5) 

Such a transformation does not change the tensor field F^ v defined in equation (3). One should note 
that the tensor field F^ v satisfies 

V T F M , + V U F T ^ + V^F UT = . (6) 

Equations (3) and (6) are similar to Maxwell equations relating the vector potential to the elec- 
tromagnetic tensor. It is therefore natural to interpret k^ and F^ v in just the same manner and 
transformations (5) as gauge transformations [19]. Hence the meson fields k^ may be recognized as 
photons. However on the basis of the axiomatic formulation of spacetime theory presented by Ehlers, 
Pirani and Schild [20], this interpretation is unacceptable. It seems that the length curvature F^ u 
must be related to a phenomena which has not been observed in nature. It has been shown that Weyl 
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mesons do not interact with leptons or quarks and other vector mesons in minimal form [21]- [23]. 
They only interact with gravitons and scalar mesons, that is, the Higgs fields. The question then 
naturally arises as to whether Weyl mesons may account for at least part of the dark matter in our 
universe. Some attempts on this subject can be found which consider the global structure of the 
universe [24]-[27]. In this paper we concentrate on the possibility of the Weyl-Dirac theory predicting 
dark matter effects on the dynamics of astronomical objects. 

3 The model 

Consider the Weyl-Dirac action functional [28] 

S[</>, k»\ = \j d 4 x^{±F^F^ + *R0 2 + a^*"} , (7) 

where a is a dimensionless constant and (j>* is the co-covariant derivative of the scalar field 4> and is 
defined as 1 

(f>*n = 4> I jl + kfj,<f> , (9) 
where *R is the modified curvature scalar given by 

*R = R + Gk^-Qk^ , (10) 

with R being the scalar curvature. The original Weyl action contains the term (*R) 2 instead of the 
term *R(j) 2 . Dirac proposed the term *R(j) 2 , thereby avoiding the great complication of the Weyl 
action. Using definitions (9) and (10) one arrives at 

S[4>, M = \j d 4 x^ I ^-F^F^ + R<j? + + (a - 6)^ 2 k^ + 2(a - 6)#^ 

+ 6(^^);4 , (11) 

where the term 6(cj) 2 k IJ ')-^^/—g is a total differential and can be ignored. 

Action (11) is invariant under local transformations of units, that is under general conformal 
transformations. Here, local unit transformations are taken to be composed of gauge transformations 
defined in equations (2) and (3) with 

9iiv — >^ 2 (av)fiV, (12) 

<\>^^r\x^ . (13) 

In order to incorporate the matter fields, we add the corresponding action S rn to (11) and consider 
S m as being built out of matter fields in the usual manner. Variation of S[cj), k^] + S m with respect 
to g^v, 4> and k^ respectively gives 

Gfiu = <j>~ 2 {T^ u + + M „ + $V> , (14) 
a<% - R<f> - (a - 6) - - ip = , (15) 

-{F^ + {a- 6)(</>V + + r = . (16) 



1 In Weyl's geometry, a localized quantity ip is called a co-tensor if under local unit transformation, it transforms as 
tp — Q," (x ^tp . Then ip is said to be of power n. Also, a co-covariant derivative is a modified covariant derivative and 
is a co-tensor. For a scalar field S it is defined as 

S*[_i — Sf_i Tikjj,S . (8) 

Here (f> is a scalar field of power n = — 1. 
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Here G^ v = R^ v — \g^ u R is the Einstein tensor and 

Tfxv = (-2 + |) g^<j) a (f> a + (2 - aWpfa - 2g^{4> a ). a + 2# w „ , 

Efu, = -F a/3 F al3 'g^ - F^ a F u/3 g a ^ , 

© M i/ = (a - 6) j</> 2 ^-k^k u + ^g^ v k a kc}j - (f> (k^ v + fc„0 M - k a <\> a g^) 
where £^ is the stress tensor of the matter fields defined as 



J /J,U 



and J 1 * denotes the current density vector 



with ip given by 



J" 



2 SSm 
/=g 5g^ ' 



5k^ ' 



5S m 



(17) 

(18) 
(19) 

(20) 

(21) 

(22) 



It should be noted that equations (14)- (16) are not independent. 

Now, taking the trace of equation (14), using definitions (17)-(19) and comparing the result with 
equations (15) and (16), one finds 



££ + (a-6) (file" + <fx^ +^ = 0, 



where from (16) we have 



(a -6) 4> 2 k» + 4>(j)» 
Taking the divergence of (14) and using (24) leads to 



+ J& = 0. 



- — S(! = J»K - \ k v + — )Jff, . 



One may assign to the non-gravitational part of the fields a stress tensor defined as 



It may be proved [30], using equation (24) that 

From equations (25) and (27) we find the following relation for the total stress tensor 



0. 



(23) 
(24) 

(25) 
(26) 
(27) 

(28) 



It is worth noting that equation (25) may be used to derive the equation of motion for a test particle 
[29, 30]. 

To progress further, let us take the matter content in our model as having consisted of identical 
particles with rest mass m and the Weyl charge q in the form of dust, so that 
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E"" = pu v u v j (29) 
where u M is the 4-velocity and the mass density p is given by 

p = mp n , (30) 

with p n being the particle density. Making use of the conservation of the number of particles, we find 
using equation (25), the equation of motion 

+ { \ ) u»u x - - = —u v F v ^ , (31) 

as 1 J m 

where j ^ j denotes the Riemann Christoffel symbol. In the left hand side of equation (31) the 
term containing ^ may be taken as a variable mass term for the particle. In such a case we have [30] 

m = mQ(j) , (32) 

where rriQ is a constant. 

4 Field of a mass source 

Let us study the field of a mass source at rest at the origin and assume that it is in a current-free 
region. In such a case we can consider the vacuum solutions of the field equations, that is we take 
^uu = and = 0. Now consider the spherically symmetric line-element 

ds 2 = -e v dt 2 + e x dr 2 + r W , (33) 

where v and A are functions of t and r. Neglecting the factor a — 6 which can be shown to be small 
[30, 10] in equation (16), we find, using the above metric 

k ,r = (34) 

where j(t) is an arbitrary function of time t resulting from the integration. Using relation (34) in 
equation (14), one finds 

e-'f-^^-^k-UtlVi, (35) 



r r 2 / r 2 (f> \ <fi 2 ) r^cf) 2 ' 



r (j) 

The form of the left hand side of equations (36) and (37) suggest that 



(37) 



e v = f(t)e~\ (38) 
Using equations (37) and (38) we see that equations (35) and (36) are identical provided f(t) satisfies 



4 = 2 ^i_ 4 ^. (39) 



The last expression gives 
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/ 



(40) 



where (3 is a constant parameter. We are now in a position to obtain the generalization of the results 
derived in [30], that is 



e = 



1 e"\ 



'2^4 



with 



—X 



Imp 2 . 
— + + A, 



A 



2 (fir <fi 2 r 2 
1 m 



+ 



2 0r 2 r 2 



+ fi3 2 (fi 2 r 2 



(41) 



(42) 



(43) 



where we have assumed ^(t) = %. The parameters m and q are constant and may be interpreted as 
mass and charge of the source. 

The above solutions are valid for any gauge function <fi(t). We choose to work in the Einstein 
gauge by taking (ds 2 — > ds 2 = cfi 2 ds 2 , — > = 1). In this case one gets 



d s 2 = - e ~ x dT 2 + e 2 ? T (e x dr 2 + rW) , 



where 



which immediately yields 



<fi = foeF. 



(44) 



(45) 



(46) 



It should be noticed that the factor <fi> 2 associated with in (14) can be considered as the gravi- 
tational coupling G in natural units (H = c = 1). Hence we may consider 

G = <fir 2 . (47) 

Using (46), the last expression leads to 

G = G e- 2fST , (48) 

where Go = <fi~ 2 may be interpreted as the Newtonian coupling constant. Therefore, the exponential 
factor e~ 2 @ T shows the time evolution of G where T is the elapsed time. From (48) we find 

t - 

The empirical measurements on the time evolution of G shows that [32] 

Gt 



G 



= (4 ±9) x 10 -13 yr -1 . 



This sets the following restriction on j3 

-6.5 x 10~ 13 < < 2.5 x 10~ 13 yr _1 . 



(49) 

(50) 
(51) 
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5 Rotational velocity of galaxies 



In this section we apply the results obtained above to the rotational velocity of typical spiral galaxies. 
In order to obtain the gravitational potential U of the point mass source m considered in section 4, 
we first note from (32) that in the Einstein gauge we have m = tuq. Now let us consider the weak 
field limit 

Soo = -(1 + 2*7). (52) 
From (44) and (42) in the conventional units, this leads to 

m G fm G\ 3 p 2 fm G\ 2 2 fm G\ 2 B 2 r 2 

where we have neglected the term \ for large r. As it is clearly seen, equation (53) differs from 
the Newtonian gravitational potential by a number of correction terms. The presence of the variable 
gravitational coupling G may equivalently account for this discrepancy, in accordance with reference 
[35]. One should note that parameter (3 appears in all the correction terms. Taking into account the 
possible bound on (3 which is given in (51), the correction terms are too small even if we consider rriQ 
to be of the order of a typical galactic mass. However, the last term, , may result in significant 
effects at galactic scales. Therefore, the first and second terms of (53) may have important effects on 
galactic and larger scales. Hence we consider the gravitational potential of a point mass source as 
follows 

TT/ . m G (3 2 r 2 

U(r) „ £_ + (54) 

This gravitational potential differs from the Newtonian one by a repulsive correction term. 

In order to study the consequences of the gravitational potential (54) on galactic objects let us 
generalize the gravitational potential U for an extended object taken as a galaxy. To this end one may 
treat the exterior of a galaxy as a vacuum solution and its interior as filled with matter and discuss 
the junction conditions along the surface where the two solutions are to be joined. Alternatively, one 
may treat the whole space as filled with matter with an appropriate fall off condition on the density. 
In fact, we adopt the second strategy in a weak field gravitational approximation. This strategy has 
also been considered in [10, 11, 12, 13, 14]. We assume that a galaxy contains galactic objects such 
as solar systems which, on galactic scales, would be seen as point like sources and the gravitational 
potential for each individual galactic object obeys the form U given by equation (54) when mo = Mq, 
where Mq denotes the solar mass. We shall consider the matter density of these objects satisfying 
an appropriate fall off condition on galactic scales. In a weak gravitational field approximation the 
gravitational potential of the galaxy as a function of r, that is, the distance from the center of the 
galaxy, is given by 



V(r) = J A7rP ^ r>) U(\r - r'\Y 2 dr> . (55) 



4Trp M (r') 

*0 

Here Pm^') is the mass density of the galaxy. We shall consider a galaxy core model for the mass 
distribution by assuming a spherically symmetric galaxy with a core density p c within a core radius 
r < r c . It should be noted that in order to more accurately describe the core behavior, one may 
consider the bulge or thin disc effects according to the observational distribution of the luminous 
matter of a typical galaxy. However, the model we consider here yields a reasonable description of 
the velocity rotation curves without taking the above mentioned effects into account. 
From (55) the gravitational potential can be obtained as follows 

GM(r) P 2 M 2 (3 2 I cm 

v{r) = — + 2Me r + m^- (56) 

In equation (56), M is the galaxy mass 
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M = M* + M HI + M DB , (57) 

with M#, Mff/ and Mdb respectively denoting the visible mass, the neutral hydrogen mass and the 
possible dark baryon mass and gas with I cm being the moment of inertia about the center of mass 
and 

M(r) = 4tt f r dr'r' 2 p c (r'), (58) 



is the mass inside the luminous core of the galaxy described by a ball of radius r = r c . Inside the 
ball the dynamics is described by Newtonian theory and the effect of the last two terms in (56) is 
negligible. Moreover outside the core the effect of the second term becomes more important. A simple 
model of M(r) is given by [14] 

M{t) = m(-^X\ (59) 



where 



j 1 HSB galaxies 

I 2 LSB and dwarf galaxies. 



Well outside the core radius, namely for r » r c , from (59) we have 

M{r) — ► M. (61) 

Also, well inside the core radius the density p(r) has a constant value for HSB galaxies and p(r) oc r 3 
for LSB and dwarf galaxies. From (56) the square of the rotational velocity of a solar object about 
the center of a galaxy can be obtained as 

2 . , M(r)G M n2 2 
1 + ^ 17 ■ ( 62 ) 

6 Observational constraint on the model parameters 

To find the consistency of this model with observational data, we compare the predictions of our 
model with data directly obtained from observations and find the best fitting parameters. Since one 
cannot expect the theory to exactly explain observational data, we give confidence intervals for the 
free parameters of the model using likelihood analysis. 

To begin with, using equation (59), we rewrite equation (62) as follows 



a. 
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v 2 ( r ) = - + a/3V . (63) 

r \r + r c J 

where a = M/M@. In order to compare the theoretical results with the observational data, we 
must compute the rotational velocities given by equation (63). For this purpose, we redefine free 
parameters of the model as 

a = 10 10 x a (64) 

= 1(T 14 x (—) (65) 



Let us compute the quality of the fitting through the least squared fitting quantity \ 2 defined by 

X \aJ,r c ) = g [Wri)-Mr i ;a > ^,r e )] 2 > (fi6) 

i % 
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where <jj is the observational uncertainty in the rotational velocity. To constrain the parameters of 
model, we use the Likelihood statistical analysis 



C(a, P, r c ) = Ne-* 2 ^~^/ 2 , (67) 

where N is a normalization factor. In the presence of every nuisance parameters, Likelihood function 
should be marginalized (integrated out). Using equation (66) we can find the best fit- values of the 
model parameters as the values that minimize x 2 (&i Pi r c)- Table I shows different priors on the model 
parameters used in the likelihood analysis. 

Table I: Priors on the parameter space, used in the likelihood analysis. 



Parameters 


prior 




a 


[0-15] 


Top-Hat 


P 


[0 - 20] 


Top-Hat 


P 


10.0 ±0.5 


Gaussian 


r c 


[0 - 7]kPc 


Top-Hat 


7 


1 or 2 


Depends to the kinds of Galaxy 



The best-fit values for the parameters of the model at 1ct(68.3%) and 2ct(95.4%) confidence inter- 
vals with corresponding reduced xl = X 2 /N {N is number of degrees of freedom) for various galaxies 
are presented in Table III. Table IV shows the best-fit values, using the Gaussian prior for (3. To 
infer the Gaussian prior, we rely on the determination of (3 from empirical measurement of the time 
evolution of G mentioned in section 4. Figure 1 compares the fitting rotational velocity curves derived 
by using the best fitting parameters of some galaxies with observational data points. 

Since a proper theory should have a constant P, we infer a mean value with its variance at la 
and 2a confidence level for (3 from the best fit-values reported in Table III and IV as 



and 




The fist value corresponds to the Top-Hat prior for (3 and the second is related to the Likelihood 
analysis with Gaussian prior. As mentioned before, due to some systematic and random errors in the 
observational data, one cannot expect the Likelihood analysis to give a unique value for (3 for our 
theoretical model so that one is inclined to find a confidence interval. On the other hand we should 
point out that in spite of many attempts in previous models [33], we have found a mostly confined 
interval for (3 which is in agrement with pervious predictions such as empirical measurement of the 
time evolution of G. In addition, Weyl-Dirac theory gives a better prediction for rotational velocity 
curves than Newtonian theory. 

Note that the plots in figure 1 give the impression as if the velocity curves tend to bend slightly 
upwards. However, this is not the case since from the gravitational potential (56) obtained for an 
extended object such as a galaxy, we may infer that there is a limiting radius for the formation 

9 



of galaxies. In fact, this gravitational potential differs from the Newtonian one by two repulsive 
correction terms. For distances near the center of a galaxy the effect of the first term is important. 
However at larger distances the first term falls off and the potential is dominated by the other two 
terms. There is a limiting distance r for which we have V(r) = and far from it the force becomes 
repulsive and the solar orbit becomes unstable. Therefore, we see that the Weyl-Dirac theory predicts 
a limiting radius for a galaxy. 

The observational Tully-Fisher relation [34] implies that v\ ut ~ GqM where v ou t is the observed 
velocity at the outermost observed radial position and M is the galaxy mass. Milgrom's phenome- 
ological law [7] or MOND model predict the Tully-Fisher relation by assuming that the mass to 
luminosity ratio, M/L, is constant across all galaxies. In Milgroms phenomenological model we have 

= GMa , (68) 

at sufficiently low acceleration a <C ao with ao ~ 10~ 8 cm/s. Relation (68) does depend on the 
magnitude of the acceleration ao, but not on the radial distance r and the rotation velocity is constant 
out to an infinite range. In contrast, our model does depend on the radius r and the time variation 
rate of the gravitational coupling. Despite the beauty of Milgrom's law in explaining flat rotation 
curves of galaxies it seems problematic to embed the theory within a comprehensive relativistic theory 
of gravity. Therefore it is not clear that the theory is as successful for explaining gravitational lensing 
of clusters and other curved space time effects. In fact strong gravitational lensing indicates a larger 
mass concentration at cluster centers than accounted for by the present form of Milgrom's theory 
[43]. 

The properties of observational rotational velocity data points for some famous galaxies used in 
our model are summarized in Table II. In order to make a comparision of our results for galaxy masses 
with other approaches we also list the mass obtained by Milgrom phenomenological MOND model 
and the Moffat relativistic MSTG model. 



Table II: 



Galaxy 


Surface 
Brightness 


M(W W M Q ) 
MOND 


M(W W M Q ) 
MSTG 


ref. 


NGC 1560 


LSB (Dwarf) 


0.59±0.05 


0.79±0.05 


[36] 


NGC 3109 


LSB (Dwarf) 


0.62±0.04 


0.78±0.04 


[37] 


NGC 55 


LSB (Dwarf) 


0.91±0.07 


1.17±0.07 


[38] 


UGC 2259 


HSB (Dwarf) 


0.55±0.02 


0.77±0.02 


[40] 


NGC 5585 


HSB (Dwarf) 


0.9±0.06 


1.17±0.07 


[42] 


NGC 247 


HSB 


1.46±0.14 


2.27±0.17 


[39] 


IC 2574 


HSB 






[41] 


Milkyway 


HSB 


10.60±0.37 


9.12±0.28 


[41] 



Column six shows the references used to extract the observational data. Column three shows 
the values M in 1O 1O M which are used to fit the velocity curves in our model. In columns four 
and five we present the galaxy masses obtained from the MOND model and the MSTG model [14]. 
Comparison of the results in column 3 with that of columns 4 and 5 shows that the masses obtained 
by the Weyl-Dirac theory are within the range obtained from the MOND and MSTG models. 

7 Inter galactic effects 

The galaxy M31 and our own galaxy (the Galaxy) are dynamically the only members of the so-called 
Local Group. Astronomical observations show that the center of M31 and the Galaxy approach each 
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other with a speed of 119 m/s [44]. A possible explanation of this phenomena was given by Kahn and 
Woltjer [45]. They showed that the effect requires the reduce mass of M31 and the Galaxy to be 100 
times larger. This discrepancy is usually accounted for by the incorporation of dark matter. Recently 
the effect has been explained by the nonsymmetric gravitational theory [12]. Here we attempt to 
explain the effect by using the Weyl-Dirac theory without resorting to the dark matter hypothesis. 

We first note that in the absence of intergalactic matter, M31 and the Galaxy form a double 
galaxy. In this case, using (56), the equation of motion for the system becomes 

,. M LG G M LG (3 2 

r H t. — r = (69) 

where Mlg is the total mass of M31 and the Galaxy and r is the relative distance between them. 
Because the distance between M31 and the Galaxy is about 700 kpc, the term M ^ G in (69) is 
negligible and can be ignored. In this case we have 

.. M LG f3 2 

rH — r = 0. (70) 

This equation is equivalent to the equation of motion derived by Kahn and Woltjer by assuming 
the intergalactic dark matter spreading homogeneously through the Local Group with mass density 
p = -^qm~P 2 - Therefore, it appears that the Weyl-Dirac theory may help us to make a better 
understanding of the intergalactic effects without the exotic dark matter hypothesis. 



8 Gravitational lensing 

The gravitational lensing or deflection of light by massive bodies provides another way to test the 
relativistic effects of gravitational theories at large scales. The gravitational lensing gives the most 
accurate information about the mass of the core of galaxies [46] and the cosmic telescope effects of 
gravitational lenses have enabled us to study faint and distant galaxies which happen to be strongly 
magnified by galaxy clusters. Also, the statistics of gravitational lensing events can be used to esti- 
mate cosmological parameters [47] . Observations of strong lensing by clusters indicate a larger mass 
concentration at cluster centers than their visible matter [2]. This discrepancy is usually accounted 
for by the existence of dark matter. Here we describe this effect by the Weyl-Dirac theory. 
The deflection angle of light rays passing through a gravitational potential V is given by 



2 f° 

= ~ / V±Vdl, (71) 

C J s 



where the integral is evaluated along the path traversed by light joining the source (s) to distant 
observer (o) and Vj_ denotes the derivative in the direction perpendicular to the path. Let us now 
consider deflection the of light from a cluster of galaxies. At these scales we can safely use the 
point mass approximation (54). Clearly, to study microlensing effects for smaller scales, we may use 
gravitational potential (56). Using (54), we find the deflection angle 

AGM f b 2 (3 2 d os \ 
a = ^ 2 T{ 1 + ^MG-)> (72) 

where b and d os denote the impact parameter and the distance between the observer and source, 
respectively. Relation (72) shows a larger light deflection than is expected from standard gravitational 
theories. For example for a cluster with b ~ 4 x 10 4 kpc , d os = 10 Gpc and (3 ~ 0.74 x 10~ 13 /yr we 
obtain 

bVd °- 10. (73) 



2M & G 

Therefore, the increase in the light deflection may be explained without the dark matter hypothesis. 
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9 Conclusions 



We have considered the Weyl-Dirac theory and showed that the spherically symmetric solutions of the 
theory lead to a gravitational potential which differs from the Newtonian one due to the appearance of 
a repulsive correction term increasing with distance. The correction term appears to be determined 
by the time variation rate of the gravitational coupling. We saw that the correction term for an 
extended massive source such as a galaxy becomes more important, especially at galactic distances. 
Hence one may expect the theory to imply significant consequences for galaxies and on galactic scales. 

We also performed well-known Likelihood analysis on the free parameters in our model to check 
the consistency of our predictions to that of the observational data. In doing so, we considered 
eight galaxies which are chosen from different types, namely the low and the high surface brightness 
galaxies. Using two different types of prior for model parameters, the best fit- values for free parameters 
at 1<t(68.3%) and 2cr(95.4%) confidence intervals have been determined. The mean value for (3 using 
Top-Hat and Gaussian priors are 6.38^34^'^ and 5.72l}'i|^2 69, respectively. However we found a 
relatively wide interval for the best fit-values of (3 for each galaxy, but this gives a small deviation in 
the goodness of the fitting curves with observation, in other words the goodness of the fit is almost 
insensitive to major fluctuations of (3. Also Likelihood analysis lets us to obtain a much more smaller 
confidence interval for (3 than the former extracted values which are also consistent with previous 
inferring. The fit to data seems to be acceptable in terms of the total mass M of each galaxy and 
a best fit-value for (3 without the exotic dark matter hypothesis at the expense of a rather wide 
interval range for (3. This could be regarded as the main disadvantage of the model presented here. 
However, it should also be pointed out that the goodness of the fit is acceptable for the underlaying 
observational data set used in this analysis. 
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Table III: Best fitting parameters at la and 2a confidence levels with corresponding 

for various galaxies. 



Galaxy 


Q 


/3 


r c kPc 






NGC 1560 


1 20+ 005 

1 20 +al0 
i - zu -0.35 


6 Q5+ 1 - 25 
°- yo -o.io 


1 oo +0 - 05 

i - uu -o.io 
1 00+ 010 

l.UU_o.20 


2 


3.26 


NGC 3109 


0-6518:* 

0.651^ 


q 40+ 100 
y.^u_ 2 10 

9 4n+ 4 - 10 

y - 4U -3.50 


1 30+ - 15 
1 S0+ ' 30 


2 


0.99 


NGC 55 


O qCr+0.30 

q qrr+0.80 


2 ™ +0 - 50 

Z - OU -0.55 

2 ^n +1 - 00 


1 60+ 005 
1 60+ - 15 

-■-•^—0.15 


2 


0.86 


UGC 2259 


o.95±8:iB 

0.95lg;» 


11.0U_ 2 30 

11 30+ 1 - 10 
ii.ou_ 460 


u - 4O -0.05 

4V KU5 
u - 40 -o.io 


2 


0.46 


NGC 5585 


2 90+ - 25 
z - yu -0.30 

2 qn +0 - 55 
z - yu -0.55 


o nc+0.70 
O.UO_ 20 

°- uo -0.80 


i 9 c+0.05 
i ' zo -0.10 

i oc;+0.10 
i - zo -0.15 


2 


2.39 


NGC 247 


9 o n +0.25 
z.ou_ 25 

9 qn+0-35 
^•<JU_o 45 


fi 35+ - 50 
6 35+ 1 - 40 

D - oo -0.90 


2.so±8:ii 
2.8oi8;i8 


1 


1.05 


IC 2574 


1 n5+ - 45 

l.UO_ 25 

1 n^ +1 - 00 


'• yo -1.50 
7 q^+3.30 


r qc-+0.65 
c qr+1.50 


1 


1.63 


Milkyway 


14.32±g:» 
14.32li;» 


q c;c+0.10 

q c;c;+0.35 
O.OO_ 25 


1 00+ 005 
1 00+ 010 


2 


1.75 
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Table IV: Best fitting parameters at la and 2a confidence levels with corresponding 
for various galaxies using the Gaussian prior for (3 in the Likelihood analysis. 



Galaxy 


a 


/3 


r c kPc 




•> 

xi 


NGC 1560 


1 20+ - 10 
1 20+ - 15 


6.95+2:12 

6.9512:1 


1 oo +0 - 05 

i - uu -0.05 
1 - uu -0.10 


2 


3.29 


NGC 3109 


1 - uu -0.10 
l.UU_ 20 


D - ' °-0.45 
D - ' °-1.05 


1 ^n +0 - 05 
i - ou -o.io 

i - ou -0.15 


2 


1.09 


NGC 55 


2 z^+ - 25 
z - 40 -0.15 

o 4 K+0.60 
z - 4O -0.35 


4 D5+ - 25 
4 - uo -0.55 

4 n5 +0 - 70 


1 40 +0 - 15 


2 


1.93 


UGC 2259 


1 - ou -0.05 
1 - ou -0.15 


°- ' u -0.50 
°- ' u -1.30 


u - oo -0.05 

u - 00 -o.io 


2 


0.79 


NGC 5585 


o i c+0.20 
z - io -0.05 

o i c;+0.45 
z - io -0.25 


4 60+ - 15 
4 60 +a70 


1 in+ - 05 

J - iu -0.05 

i.iol2:i2 


2 


3.29 


NGC 247 


2 3f)+ - 15 

9 on+0.35 
^- OU -0.35 


fi 35+ - 45 
fi S5+ 1 - 05 


2.8oi2i2 
2.8oi2:i 


1 


1.05 


IC 2574 


1 40+ - 20 


6-6012:1 

6.601^ 


r nc;+0.70 
°. yo -0.65 

5.95l2: 6 7 2 


1 


1.65 


Milkyway 


i3.4stS:?i 

13.45±J:I° 


3-80l2i2 
3.80l2i2 


o.95l2:l 
o.95l2:12 


2 


2.09 
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NGC274 NGC 55 




Milkyway UGS 2259 




Figure 1: The rotational velocity of various galaxies in terms of distance. The solid line corresponds to the theoretical 
prediction, equation (63), with the best fit values for the free parameters derived by using Top-Hat prior for fl. The 
theoretical fitting corresponding to the Gaussian prior for (3 has been indicated by the dashed-dot curves. The dashed 
curve shows the Newtonian-Keplerian prediction. Data from observations indicated by symbols are for the following 
galaxies: NGC 247, NGC 55, NGC 1560, NGC 3109, IC 2574, NGC 5585, UGC 2259 and Milkyway. 
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